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ABSTRACT 



Aims. We investigate the photometric signature of magnetic flux tubes in the solar photosphere. 
Methods. We developed two dimensional, static numerical models of isolated and clustered magnetic flux tubes. We investigated 
' the emergent intensity profiles at different lines-of-sight for various spatial resolutions and opacity models. 
Results. We found that both geometric and photometric properties of bright magnetic features are determined not only by 
the physical properties of the tube and its surroundings, but also by the particularities of the observations, including the 
line/continuum formation height, the spatial resolution and the image analyses techniques applied. We show that some 
observational results presented in the literature can be interpreted by considering bright magnetic features to be clusters of 
smaller elements, rather than a monolithic flux tube. 
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1. Introduction 

Recent observations taken in the G-band with sub- 
, arcsecond resoluti on have al lo wed the study o f pho- 
tometric feerger et al. I 2004 : Okunev fc Kneer I2OO4I : 



namic 



Ishikawa e^al. I l2007l: Jgerger et al.l l2007t 



Rouppe van der Voort et al. 1 120051 : llshikaw a et al. 



while Hirzberger fc Wiehr ( 20051 ), after analysis of 
extreme limb observations, reported no trend in contrast 
or size with position on the solar disk. Discrepant CLVs 
have also been obtained when compa ring G-band with 



iBovelet fc Wiehr I 1200^; iNisenson et al. 




other continua measurements 



1991; Sutterlin 




(e.g. 



Auffret fc Muller 



Sanchez Cuberes et al. I 12002: 



properties of small magnetic flux concentrations in the so- 
lar photosphere. Observations away from disk center have 
revealed that bright magnetic flux concentrations are usu- 
ally accompanied by a dark lane on their dis k center side 
and by a b right ta il on the limb- ward side (iLites et al~l 



I2OO4I : IHirzberger fc Wiehrl l2005l : iBerger et al.l l2007h . 



However, some bright points show no associated tail or 
dark lane, and some dark lanes with enhanced magnetic 
flux d ensity show no associated bright point (Berger et al. 
|2007[ ). Somewhat contradictory measurements have been 
reported for the center-to-limb variation (CLV) of G-band 
bright point p hotom etric contrast and size. For example 
iBerger et al.l ( 2007t ) reported maximum contrast at a 
position much closer to the limb than previous observers. 



Okunev fc Kneer I |2004| ). Together these observations 
have raised questions about the validity of viewing 
bright magnetic features as monolithic flux elements. 
In particular, it has been suggested that the observed 
CLVs of contrast and size are determined by properties 
of u nresolved clusters of magnetic elements at the limb 
(e.g. lOkunev fc Kneer II2005I : IBerger et al.ll2007[ ). 

The properties of small individual magnetic fea- 
tures have been intensively investigated using both 
two and three dimensional magneto-hydrodynamic 
numerical simulati o ns (e.g. iKeller et al. I 2004 : 
'Carls son et aTl l2004l: Isteiner M2OO5I : iHasan et alH l2005l : 
Uiten broek fc Tritschlerl l2006h . Some of the highlights 
of these modeling efforts include the following. Using 
both st at ic an d dynamic two-dimensional simulations, 
( 2OO5I ) showed that the dark lane observed 
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inside and outside) the 
magnetic tube. Simil ar 



disk-center-ward flank of the 
results were previous l y pre - 



sented by [Pizzo et al. I (Il993al lbl): IPeinzer et aTl (|l984l ): 



Knolker et al. I ('l988a''b)~ who also investigated the 
sensitivity of the contrast intensity profiles to the physical 
properties of the tube (such as the temperature boundary 
condition, the evacuation of the tube, the intensity of 
magnetic field, and the inhibition of convection) and its 
temporal evolution. Ultimately, the intensity contrast 
profiles are determined by the degree of evacuation of 
the tube and the properties (particularly temperature) 
of the atmosphere inside and surrounding it. Fully 
three-dimensional ma^ne t o-hydrodynamic simulations 
(e.g. Keller et al.1 l2QQ4l: ICarlsson et aTl 2004) have 
demonstrated that the increase in magnetic element 
contrast with distance from disk center largely reflects 
the properties of the solar granulation behind the tube. 
This becomes visible away from disk center because of 
the reduced opacity within t he tube for highly oblique 
viewing angles. Moreover, iDe Pontieu et al. ] 12006), 
using both magneto-hydro-dynamic simulations and 
high spatial resolution observations, demonstrated that 
the temporal evolution of the photometric properties of 
magnetic flux tubes occurs on the characteristic granular 
time scale. 

The literature addressing the properties of fiux tube 
aggregations is less extensive. iCaccin fc Severino I (|l979h 
showed that the shape of the contrast- CLV of clus- 
ters of flux tubes is independent of the shape of 
the single underlying model tube (fu n nels or cylin- 



der s)^^lore_j;ecentlx2JP^^ fc Kneer 
ing iKarpinskv fc Okune v (1998), used three-dimensional 



( 2005h . follow- 



static models in pressure equilibrium, to show that the ap- 
pearance of aggregations of flux tubes and their measured 
polarimetric signals are strongly influenced by their posi- 
tion on the solar disk, filling factor and spatial resolution. 

In this paper we investigate how the photometric sig- 
natures of isolated and clusters of small magnetic elements 
(such as the center-to-limb variation of photometric and 
geometric properties) reflect, not only the physical prop- 
erties of the tube and its surroundings (such as temper- 
ature stratification and magnetic flux density), but also 
the particularities of the observation itself (such as spatial 
resolution, formation height of the wavelength of observa- 
tion and image analyses technique employed). We examine 
these sensitivities using two-dimensional models of static 
magnetic flux tubes in radiative equilibrium with their 
surroundings. We show that results obtained by the anal- 
yses of recent high resolution observations can be inter- 
preted assuming that observed magnetic bright features 
are actually aggregations of unresolved magnetic struc- 
tures. 

In §2 we describe the model. In §3 we investigate ther- 
mal properties of aggregations of magnetic flux tubes. In 
§4 we present the obtained contrast profiles and discuss 
their dependence on observational wavelength and spatial 
resolution. In ^5 and ^6 we show how the observed CLVs of 



magnetic features depend not only on peculiarities of ob- 
servations, but also on the numerical techniques employed 
to detect features on images. In §7 we investigate the size 
dependence of contrast. Conclusions are presented in §8. 

2. The Model 

We consider static models of isolated and clustered mag- 
netic flux tubes in Radiative Equilibrium (RE, hereafter) 
with the surrounding non-magnetic plane parallel atmo- 
sphere. The RE is imposed with an iterative sch eme sim- 
ilar to the one proposed bv IPizzo et al. I ( 1993ah In brief, 
initial atmospheric conditions for the quiet sun and flux 
tube are imposed, the mean intensity everywhere in the 
domain is estimated and the logarithmic gradient V = 
dlnT/dlnP, where T is the temperature and P is the 
pressure, is calculated. A new temperature field is then 
computed imposing RE in those layers which do not satisfy 
the Schwarzschild criterion. A smooth variation of tem- 
perature between convectively stable and unstable layers 
is imposed using a linear interpolation. Namely, in the 
region 70 km thick where V ^ V ad^ the temperature is 
given by T(x, z) = (I- a)T{x, z)r + aT(x, z)c, where Tr 
is the temperature value estimated by imposing RE, Tc is 
the temperature value given by the assumed atmosphere 
model and a is a parameter which varies between and 1. 
From the solution of the hydro-static and state (a perfect 
gas is assumed) equations the new pressure and density are 
computed respectively, while source function is estimated 
assuming local thermodynamic equilibrium (LTE). In this 
new atmosphere the mean intensity and the logarithmic 
gradient V are re-evaluated. The procedure is halted when 
temperature relative difference among two consecutive it- 
erations is less than 0.001%. The resulting atmosphere 
therefore differs from the initial one only in those layers 
which are stable against convection. The mean intensity 
i s evaluated in tegrating over the solid angle by a Carlson 
( Carlson I Il963 ) quadrature scheme (namely, scheme A 
with 10 angles per octant). The intensity in each spatial di- 
rection is evaluated by a code (ICriscuoini2007lil| ba sed on 
the short-characteristic technique (jKunasz fc Aiier ][l988,) . 
The initi al non magnet ic atmosphere is derived from the 
model of lKuruczl (|l994 ): the opacity is assumed to follow 
a power law in pressure and temperature . For a complete 



description of the atmosphere model see iGiordano et al. 
(|2003). 



The presence of a flux tube is approximated by shift- 
ing a portion of the atmosphere downward prescribing the 
initial Wilson depression amplitude, which along with the 
size of the tube is a free parameter of the simulations. 
Note that, when iteration is halted, the final Wilson de- 
pression differs from the initial one. Incoming radiation at 
the top of the two-dimensional domain is set to zero while 
radiation at the bottom of the domain is estimated based 
on the radiative diffusion approximation. Horizontal peri- 
odicity is imposed at the vertical boundaries. The spatial 
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domain is sampled by a grid with vertical and horizontal 
resolutions of 7.1 km and 3.5 km, respectively. The verti- 
cal domain is 1130 km deep, 630 km above and 500 km 
below the optical depth unity. The horizontal extent of the 
domain is such that the ratio between the tube diameter 
(or of the cluster) and the width of the domain is large 
enough so that the presence of the tube does not signifi- 
cantly influe nce the temperature stratificati on of the quiet 
atmosphere ( Fabiani-Bendicho et al. Ill992 ). 

The model solutions we present are thus based on the 
following assumptions. The geometry we consider is 2D, 
although the radiative field and quadrature employed is 
fully 3D. We expect this to reduce the effects of radiation 
channelling (see next section), although temperature dif- 
ferences between 2 and 3p calculations are estimated to 
be small ('Uiten broek lllQQSl ). Our model is also not in hori- 
zontal pressure equilibrium, which means that the tube ge- 
ometry is simplified; flaring with height is not reproduced. 
As a consequence, the tube contrasts near the limb are 
somewhat reduced. However, thin flux tube models sug- 
gest that this flaring is negligible at photospheric heights 
we cons i der. F inally, our model assumes a gray opacity. 
Vogler I ( 2004 ) showed that this assumption, by neglect- 
ing addition al coupling of the radiat ion with matter (see 



for instance Isteiner fc Stenfio1ll990l ). leads to an under- 



estimation of the radiative illumination effect (see next 
section) in flux tubes, and therefore to under-estimation 
of temperature in magnetic features. That author also 
showed that magnetic features contrast is under-estimated 
in gray calculations. We expect the same effects in results 
obtained with our calculations. 

The models we developed are not meant to exactly 
reproduce the results obtained by recent high resolution 
observations, but rather to show qualitatively how mea- 
surements of contrast and geometric properties of isolated 
and clusters of magnetic elements are affected by observa- 
tional issues. 

3. Thermal properties of clusters of flux tubes 

In this section we discuss the thermal properties of clusters 
of magnetic fiux tubes, solutions in which more than one 
tube is in RE with the surrounding atmosphere. While we 
have investigated the thermal signatures of a range of flux 
tube sizes and configurations (i.e. number, strength, and 
separations), we present here only the results for aggrega- 
tions of flux tubes with diameter 70 km, center separation 
of 105 km, and initial Wilson depression 150 km. 

Examples of the temperature fields obtained by sim- 
ulations are illustrated in Fig. [TJ This shows the tem- 
perature field for an isolated (top), two (middle) and 
three (bottom) magnetic flux tubes. The final Wilson de- 
pression of the tubes in RE is approximately 60 km in 
all the models. Using the thin flux tube approximation, 
the magnetic field intensity along the tubes axis esti- 
mated at optical depth r = 1 in the non magnetic at- 
mosphere is approximately 1.5 kG in all the configura- 
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Fig. 1. Detail of the temperature field inside and around 
an isolated (top), two (center), three (bottom) magnetic 
fiux tubes in RE. Diameters are 70 km and tubes axis are 
105 km apart. The continuous lines represent temperature 
iso-contours (in Kelvin). The dotted lines represent the 
heights at which r = 1 for various lines of sight: /i = 0.2, 
0.4, 0.6, 0.7, 0.9, 1 from top to bottom. The colored lines 
in the bottom panel represent the r = 1 depth at /i = 1 
for various opacity models. Red dotted line: gray model. 
Blue continuous line: A = 8000 A. Green continuous line: 
A = 5000 A. Red continuous line: A = 16000 A. 



main is lower then the surrounding atmosphere, as im- 
posed by b oundary and initial conditions. Because of the 
channeling (Cannon Ill970l: Kneer fc Truiillo-Bueno |[l987l : 



iFabiani-Bendicho et al. 1 19921 : Hasan et al. 1 19991 ) of radi- 
ation escaping from the surrounding hotter atmosphere, 
the temperature difference between the magnetic elements 
and the non magnetic atmosphere decreases with decreas- 
ing optical depth. Eventually, in the optically thin part 
of the domain the temperature of the tubes and of an 
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decrease o f tempe r ature in clusters of tubes was predicted 



by Dein zer et al. I (19841). 



Fig. 2. Temperature along central axis versus optical 
depth for isolated, three and five tubes configurations in 
previous figures. The temperature profile of non magnetic 
atmosphere is given for reference. Arrow indicates optical 
depth unity along the axis of the isolated tube model. 



of the non magnetic atmosphere. In particular, we note 
that between 0.5 < r < 2 the temperature of plasma 
surrounding the tubes is lower with respect to the quiet 
atmosphere, while it increases at shallower depths. The 
depth at which the channeling is effective in heating the 
tubes, as well as the area of the surrounding atmosphere 
affected by the presence of the tubes, is determined by 
the ratio of the tube diameter and of the horizonta l opti- 
cal depth value (jPizzo et al. lll993bl : ICriscuoh ll2QQ7h . The 
heating is larger along tubes flanks, and the effect is larger 
for external flanks of peripheral tubes (see for instance the 
isotherm at 6300°K in the three tubes case of Fig. [3]). 

Inspection of temperature fields also reveals that the 
amount of heating and cooling is a function of the num- 
ber of fiux tubes. This is better illustrated by the plot 
in Fig. O which shows the temperature variation with 
depth along the axis of an isolated tube and of the central 
tube of clusters of three and five elements. The temper- 
ature of the quiet atmosphere is given for reference. At 
heights 1 < log(r) < 2 heating is larger for aggregations 
of tubes with respect to the isolated one, but the increase 
of temperature is almost independent of the number of 
tubes in the cluster (three or five in this case). This is 
due to the fact that at these heights, because of the high 
value of the opacity (or of the horizontal optical depth), 
radiation escaping from the furthest tubes cannot pene- 
trate into the central tube and contribute to its heating. 
On the contrary, at r < 1 heating is largely dependent 
on the number of tubes. In these regions the increase of 
temperature is determined by the " illumination" from the 
layers at r = 1 wit f iin and along the flan ks of the tube 



( Pizzo et al. Ill993bl : iKnolker et al. Ill988a[ ). The increase 



of the number of magnetic elements increases the width of 
the illuminating surface and therefore increases the tem- 
perature in the higher layers of the photosphere above the 
cluster. Between < log(r) < 1 temperature is lower in 
clusters than near isolated tubes. There is a little depen- 



Finally, we investigated the sensitivity of the temper- 
ature field to the increase of the number of adjacent flux 
tubes for various magnetic field intensities (for various val- 
ues of the initial Wilson depression). We found that, when 
holding the number of tubes constant, the stronger the 
magnetic field, the greater the increase of temperature at 
r < 1 and the larger is the temperature decrease between 
< log(r) < 1. 

These results a re in agreement with t h e one s previ- 



ously obtained bv iFabiani-Bendicho et al. I (|l992l ). Those 



authors investigated the dependence of temperature vari- 
ations on filling factor by numerical 2D simulations, ex- 
ploiting the horizontal periodicity of the intensity field 
and varying the horizontal size of the spatial domain of 
the simulations. Their models show the same temperature 
increase in the outer layers and decrease in lower layers of 
the atmosphere. 

4. Contrast profiles 

We define the contrast to be the ratio I{x)/Iq — 1, where 
I{x) is the emergent intensity value at each horizontal po- 
sition X and Iq is the emergent intensity value far from the 
tube. We find, for various flux tube models and clusters 
configurations, that the contrast profiles depend on both 
the number of tubes and their separation distances. For 
example. Fig. [3] shows the contrast CLVs of an isolated 
flux tube and groups of two, three and five flux tubes 
of 70 km diameter and separation distance of 105 km. 
At disk center, the profiles are characterized by two pos- 
itive contrast "humps" and a surrounding negative con- 
trast "ring". The "humps" are generated by the differ- 
ent inclinations of the r = 1 contour with respect to the 
isotherms as shown in Fig. [T] (see also lPeinzer et al. Ill984t 
Knolker et arill988al lbl). The "dark ring" is generated by 
the cool material surrounding the tube. We note that the 
maximum contrast value at /i = 1 slightly increases with 
the number of tubes. This reflects the slight increase of 
temperature with the number of adjacent tubes shown in 
Fig. [2] at the r = 1 depth inside the tubes. Figure [3] also 
shows that the tubes are fully resolved only at disk cen- 
ter. The shallower the line of sight, the more the tubes 
appear as a single structure, with dimming on the disk 
center side and brightening on the limb side, as observed 
for singl e fiux tubes. Thi s is in agreement with the find- 
ings of lOkunev fc Kneer (|2005l ). who showed that at the 



limb aggregations of simulated faculae appear as chains of 
small bright points. For lines of sight at which the tubes 
are still distinguishable, the ones which are closer to the 
limb appear brighter than the tubes closer to the disk. 
This is due to the fact that the radiation emitted at the 
limb side crosses more than one tube and is therefore less 
attenuated, or, in other words, the r = 1 contour pene- 
trates deeper into the atmosphere, thus sampling higher 
temperature layers, as also shown in Fig. [H For similar 
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Fig. 3. Contrast profile across an isolated, three and five tubes observed at different positions on the solar disk. Tubes 
are 70 km wide and separated (measured from tubes axis) by 105 km. The disk center is on the left. Continuous line: 
/i = 1, dotted line: /i = 0.9; dashed line: ja = 0.7; dot-dashed line: /i = 0.6; dot-dot-dashed line: fi = 0.4 
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Fig. 4. Left: contrast profile across two tubes 70 km wide 105 km apart. Right: contrast profiles of an isolated flux 
tube 175 km wide. Legend is as in figlH 



crease with the number of tubes in the cluster. Since for a 
given inclination and cluster configuration there is a finite 
number of tubes crossed by each ray, contrast increases 
with the number of tubes until this maximum number is 
reached. This maximum number increases with the incli- 
nation of the line of sight. 

Figure [4] displays solutions for isolated and double flux 
tubes which occupy the same physical space: 175 km. At 
disk center the contrast of the isolated large flux tube is 
lower with respect to the two small tubes case, while the 
contrary is observed off disk center. The opacity of the 
larger single flux tube is smaller than that of the two small 
tubes, and the r = 1 contour penetrates deeper into the 
hotter non magnetic atmosphere. The consequent contrast 
profiles at different positions on the solar disk are quite 
different for the two configurations up to approximately 
/i = 0.6. In particular, we note that the disk-center- ward 
dark lane is deeper and larger for the larger isolated flux 
tube than it is for the cluster. At smaller values of /i the 
contrast profiles of the two configurations are quite simi- 
lar. 

4.1. Dependence on opacity model and spatial 
resolution 

We have computed the contrast profiles of isolated and 
clusters of flux tubes after varying the model opacity. This 
has allowed us to investigate the dependence of observed 



tions. Continua opacities in three different wavelengths, 
5000, 8000 and 16000 A, were computed by using classical 
formula for continua processes and scattering of princi- 
pal elements (see for instance Stix 2002). No blanketing 
line effects were considered. Abundances were taken from 
Grevesse et aTl (|l993[ ). 

We find that, because of the heating of the upper layers 
of the flux tubes, the contrast increases with the height of 
the r = 1 depth at the wavelength being considered. The 
increase of contrast is larger at disk center, where larger 
variations of the Wilson depression occur, as shown for 
instance for the cluster of three tubes in Fig. [H Figure [5] 
shows the corresponding contrast profiles at various lines 
of sight. We note that the negative contrast area surround- 
ing the tubes at disk center, which is present for the cor- 
responding gray atmosphere model illustrated in Fig. [3l 
is significantly reduced at the three wavelengths investi- 
gated. In fact, as shown in FigiU in the grey model the 
r = 1 line forms at a height at which the temperature of 
the plasma surrounding the tubes is lower with respect to 
the one of the non magnetic atmosphere. At wavelengths 
which form at larger and shallower depths the temper- 
ature differences between the area surrounding the tube 
and the quiet atmosphere are reduced and the dark ring 
is not observed. Variation with wavelength largely affects 
the contrast values for vertical lines of sight, as can be seen 
by comparison of contrast profiles illustrated in Fig. [5l In 
particular, the absolute contrast of dark lanes increases 
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Fig. 5. Contrast profile across three tubes 70 km wide 
105 km apart for various lines of sight and three different 
wavelengths. Blue, dot-dashed line: A = 8000 A. Green, 
continuous line: A = 5000 A. Red, dashed line: A = 16000 
A. The disk center is on the left. From right to left: /i = 1, 
/i 0.9 and /i 0.4. 



optical grey depths larger than one), while they gradually 
disappear at wavelengths which sample shallower layers of 
the atmosphere (i.e. at optical grey depths smaller than 
one). Finally, the size of the bright area of the profiles (see 
next section) observed off disk center increases with the 
height of unity optical depth. 

In order to investigate the effects of spatial resolution 
on observations, simulated contrast profiles have been pro- 
jected onto the plane of the sky, that is the contrast pro- 
files have been re-sized by the factor /i, and convolved with 
Gaussian functions of different Full Width Half Maxima. 
We found that projection reduces the maximum contrast 
and the amplitude of the limb ward tail, thus producing 
smaller and more symmetric profiles. This is illustrated 
for instance by comparison of plots in Fig. [6] with those 
of Fig. [3l These show the contrast profiles of various flux 
tubes configurations after and before projection, respec- 
tively. Reduction in spatial resolution causes decrease of 
maximum contrast, broadening of profiles and smearing of 
small scale features, such as the dark lanes and the dou- 
ble "humps". These effects are larger for tubes observed 
at vertical lines of sight, and smaller for shallower lines 
of sight, where the intrinsic profiles are already blended 
by line of sight effects (Fig. [6]). In particular, due to the 
presence of dark features, at vertical lines of sight the de- 
crease of contrast values due to the reduction of spatial 
resolution is larger than at shallower lines of sights. As 
will be discussed in next section, these variations change 
the shape of the CLV. It is also interesting to note that in 
the models shown in Fig. [6] with resolution 0.1", the dark 
lane is present in the case of the isolated large fiux tube, 
but is not present for clusters of small magnetic elements. 
Inspection of Fig. [6] reveals that the effects of reduction of 
resolution also depend on the number of flux tubes in a 
cluster. For instance, maximum contrast observed at disk 
center is approximately the same for the isolated, two and 
three tubes clusters at the resolution of the simulations. 
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Fig. 7. CLV of maximum contrast for various flux tubes 
models and clusters at resolution 0.1". Continuous: sin- 
gle flux tube 70 km wide. Dotted: two flux tubes 70 km 
wide 105 km apart. Dashed: three flux tubes 70 km wide 
105 km apart. Dot-dashed: five flux tubes 70 km wide 
105 km apart. Dot-dot-dashed: single flux tube 175 km 
wide. Dashed-thick: single flux tube 280 km wide. Initial 
Wilson depression is 150 km for all the models. Continuous 
red thick line: CLV of average maximum contrasts of all 
flux tubes configurations computed imposing a contrast 
threshold of 0.08. 




Fig. 8. CLV of maximum (black) and mean (red) con- 
trast of a fluxtube 175 km wide. Spatial resolution is 0.1" 
(continuous), 0.3" (dotted) and 0.6" (dashed). Circles and 
squares indicate the maxima of the curves. 



trast proportional to the flux tube number is observed for 
the three configurations. 



5. CLV of maximum and mean contrast 

We have estimated the CLV of maximum and mean con- 
trasts of various fiux tubes models and configurations. The 
mean contrast is defined as the average of contrast values 
larger than a specified threshold. Since for aggregations 
of flux tubes several distinct features can be selected in 
this way, we considered the average value of the contrast 
of each feature. Variations of the CLVs of mean contrast 
with variation of the threshold value were also investi- 
gated. Since the application of an intensity or contrast 
threshold is usually employed to detect features on im- 
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Fig. 6. Contrast profiles across flux tubes projected onto the plane of the sky for different resolutions and positions 
on the solar disk. Top Left: single flux tube 175 km wide. Top right: single flux tube 70 km wide. Bottom left: two flux 
tubes 70 km wide 105 km apart. Bottom right: three flux tubes 70 km wide 105 km apart. Continuous line: resolution 
of simulations. Dashed line: resolution 0.1". Dot-dashed line: resolution 0.3". In each panel, from right to left: /i = 1, 
/i = 0.9, /i = 0.6, and /i = 0.4. 
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Fig. 9. CLV of maximum contrast of a cluster of three 
flux tubes 70 km wide 105 km apart observed with spatial 
resolution 0.1" in various continua. Continuous line: grey 
atmosphere. Dashed line: A = 8000 A. Dotted line: A = 
5000 A. Dot-dashed line: A = 16000 A. 

the measured contrast CLV on the identification method 
adopted. 

We found that the CLV varies significantly with the 
number of flux tubes in a cluster. This is illustrated by 
Fig.[7j in which the CLV of maximum contrast, in the case 
of spatial resolution 0.1", is plotted for various models of 
isolated and clustered flux tubes. The plot shows that the 
curves flatten and the peak moves toward the limb as the 
number of tubes increases. The shift of the peak with an 
in crease in filling fact or was also found in 3-D simulations 
bv lOkunev fc Kneer I ([2005). In particular, it is interesting 
to note that the curve in Fig. [3 obtained for a cluster of five 



CLV obtained by those authors for a cluster of elements of 
radius 100 km and filling factor 0.2. The plot shows clearly 
that contrast of large monolithic structures (isolated flux 
tubes) is lower than the contrast of smaller structures at 
disk center, while the opposite is observed off disk center. 
A sharp decrease from the center to the limb is observed 
for the isolated smallest structures, while smoother CLVs 
are observed for clusters of small elements. This flnding 
indicates that at the limb clusters and large tubes are 
more likely to be observed rather than isolated small tubes 
and therefore weigh more in the estimation of the CLV of 
contrast derived by observations. 

To investigate this effect we computed the average con- 
trast at each position on the solar disk of the various mod- 
els shown in the flgure. To mimic selection effects asso- 
ciated with the application of intensity threshold criteria 
on images, average has been computed discarding contrast 
values smaller than 0.08. The obtained CLV is represented 
by the red line in the plot of Fig. [3 The resulting CLV 
is flatt er, espe c ially a t the limb, than the CLV of isolated 
tubes. Spruit (|l976l ) investigated the CLV of maximum 
contrast obtained connecting the maxima of the CLVs of 
various tubes models. The author considered four classes 
of models and, in agreement with our flndings, showed 
that for three classes of models the curves so obtained 
were flatter than curves of individual tubes models, al- 
though the shapes of the curves were very different from 
the ones we have obtained. 

As already noted for the contrast proflles, we found 
that the maxima of both maximum and average contrasts 
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Fig. 10. CLV of size of different flux tubes models and 
clusters for contrast tliresliold value 0.08. Legend as in 
Fig. El 




Fig. 11. CLV of size of a tube 175 km wide for various 
spatial resolutions. Square: 0.1". Triangle: 0.3". Diamond: 
0.6". Contrast threshold value is 0.02. 



6. Size and asymmetry 



tion for models in which dark lanes are large and dim. This 
effect is reduced in the case of aggregations of small flux 
tubes, for which dark lanes are less pronounced. Figure 
[8] shows for instance the CLV of maximum (black lines) 
and mean (red lines) contrast in the case of an isolated 
tube 175 km wide with initial Wilson depression 150 km. 
Different line-styles represent different spatial resolutions. 
Squares and circles indicate the maxima of the curves, 
which shift toward the limb at the decrease of resolution. 
Here we define as mean contrast the average of contrast 
values larger than 0.02. The plot also shows that, for a 
given spatial resolution, maximum and mean contrast val- 
ues are correlated, i.e. the shapes of the curves are similar 
although the contrast values are different. We also notice 
that with the decrease of resolution the curves flatten and 
the decrease of contrast toward the limb is less abrupt. 
The dependence of the cont rast- CLV on spatia l resolu tion 
was previously reported bv lOkunev fc Kneer I ( 2005h . al- 
though a systematic shift of the peak with the decrease 
of resolution was not observed. This is most likely due 
to the lack of negative contrast area (dark lanes) in their 
simulations. 

We also investigated the dependence of CLV of the 
mean contrast on the contrast threshold value applied. 
We found a change in the contrast values measured, but 
not substantial variation of the CLV shapes. Finally, we 
found that the contrast- CLV is significantly dependent 
also on the opacity model, that is on the wavelength of 
observations. Figure [9] plots the CLV of maximum con- 
trast of a cluster of three flux tubes 70 km wide and 105 
km apart for the opacity models previously discussed. The 
largest changes in contrast occur at disk center (contrast 
decreases at wavelengths which form in the deepest lay- 
ers), while at the extreme limb, contrast values are nearly 
independent of observation wavelength. As a consequence, 
the position of the peak of the CLV shifts from /i ^ 1 at 
wavelengths which form in the highest layers to lower val- 



We have investigated the geometric properties of isolated 
and clustered flux tubes. This analysis was performed on 
contrast profiles projected onto the plane of the sky and 
convolved with gaussian functions of different widths as 
explained in the previous paragraph. 

We define the size of magnetic features as the largest 
distance between the points over which the contrast ex- 
ceeds a given threshold. We find that the shape of the 
size- CLV is very dependent on the flux tube models in- 
vestigated, as illustrated in Fig. [TOl In general, CLVs of 
isolated or clusters of tubes which present marked double 
humps and dark lanes, show a minimum at approximately 
/i = 0.9 and a peak between ja = 0.6 and ja = 0.4. An ex- 
ample is given by dot-dot-dashed line in Fig. [TOl which 
represents the size- CLV of a tube 175 km wide and ini- 
tial Wilson depression 150 km. The decrease of size from 
disk center to /i 0.9 is given by the smooth disappear- 
ance of the disk-centerward side wall and by the appear- 
ance of the dark lane, whose presence reduces the size of 
the "bright" area. The rise of size at intermediate lines 
of sights is due to the increase in width of that portion of 
the r = 1 contour, which penetrates into the non magnetic 
atmosphere, as shown for instance in Fig. [3l The steep de- 
crease at the extreme limb is due to both foreshortening, 
and a decrease in the size of the region with contrast val- 
ues larger th an the th r eshold applied. A CLV of this type 
was found bv lSteiner I (|2005h for a static flux tube 170 km 
wide and Wilson depression 150 km. For models which do 
not present marked dark lanes or humps at disk center we 
found monotonic variations of size with position on the 
solar disk, as illustrated by most of the curves in Fig. [TOl 
The other exception in that figure is the CLV of the 280 
km wide tube, for which an increase from the center up to 
/i = 0.4 is observed. This is due to the fact that the con- 
trast profile at disk center of this model is characterized 
by two positive contrast humps separated by a dark lane. 
According to our definition, the two humps are detected 
as two distinct features whose average size is the estimate 
of the size of the tube. This definition underestimates the 
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Fig. 12. CLV of size of a tube 70 kni wide at resolu- 
tion 0.1" for different contrst threshold values. Diamonds: 
0.005. Square: 0.02. Triangle: 0.04. Asterisks:0.06. 
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Fig. 13. CLV of skewness at various spatial resolutions 
for a flux tube 175 km wide. Diamonds: original. Plus: r 
= 0.1". Asterisks: r = 0.3". Triangles: r = 0.6". 



tube 175 km wide. We found that, in general, the skewness 
value is quite low, thus indicating that the contrast profiles 
are, overall, quite symmetric. The observed asymmetries 
are larger close to disk center and decrease toward the limb 
for the best spatial resolution investigated. Nevertheless, 
the decrease of resolution mostly affects vertical lines of 
sight, therefore for intermediate resolutions the asymme- 
tries have a peak at approximately 0.4 < ji < 0.6. At the 
lowest resolution investigated the asymmetries decrease 
again monotonically from disk center to the limb. It is 
also worth noticing that at /i < 0.4 profiles are nearly sym- 
metric. This is not in disagreement with findings of pre- 
vious authors, who reported an increase of the limb ward 
tail with the inclination of the line of si^ht ( e.g SteineT] 
I2OO5I : IPizzo et al. Ill993bl : iDeinzer et al. Ill984l ). since here 
we are considering profiles project onto the plane of the 
sky and, as previously shown, projection mostly affects 
the extreme l imb. Moreover, this is in agreement with 
iBerger et al.l ( 20071) , who found that experimental aver- 
age contrast profiles at /i = 0.4 are more symmetric than 
average profiles at ji = 0.6. We also found that the value 
of the asymmetries defined in this way increases with the 
size of isolated flux tubes and with the number of tubes 
in a cluster. Finally, it is also important to notice that 
the skewness poorly describes the asymmetries of distri- 
butions which have multiple peaks and is therefore not a 
good indicator of the asymmetries for profiles of isolated 
large flux tubes at /i > 0.9, or in the cases for which the 
profiles of clusters of tubes show multi-peaks. 

7. Contrast-size relation 



of the CLV illustrated. Finally, the red line is the average 
value of the sizes obtained by the models in the figure for 
each line of sight investigated. As for the contrast, the av- 
erage size-CLV is flatter than the single models, especially 
at the limb. 

The size-CLV curves also assume quite different shapes 
with variation in resolution. In general, the peak of the 
curve shifts toward disk center and the decrease of size 
toward the limb is steeper at lower spatial resolutions (Fig. 
[TT]) . Moreover, while a decrease in the intensity threshold 
used to define structures corresponds an increase of the 
measured size, the CLV profiles remain similar to the ones 
obtained at lower resolution. An example of these effects 
is given in Fig. [121 which shows the size-CLV of a tube 70 
km wide for four values of threshold contrast. It is worth 
noting that for high contrast threshold values or at poor 
spatial resolutions, small faint features are not detected. 

As pointed out previously, projection onto the plane 
of the sky and reduction of resolution modify the shape of 
contrast profiles. To quantify these variations we analyzed 
the skewness of the profiles and defined as "asymmetry" 
the skewness of distribution of contrast values larger than 
a certain threshold. An example of the CLV of the asym- 
metries for threshold contrast value 0.001 and various spa- 



Scatter plots of measured contrast against mag- 
netic feature size presented in the literature usu- 
ally show a large dispers ion. Fits to these points 



have s hown no trend (IBerger et al. 1 19951 : IWiehr et al. 



l2004t ISanchez Almeida et al. 
ffiovelet & Wie^ 



crease 



- i2008h 
l2Q03l 



or a 



slig:ht 



m- 



Hirzberger fc Wiehr 



2OO5I : IBerger et al.l 12007) of contrast with magnetic fea- 



ture size. In particular, IBerger et al.l (|2007l ) showed that 
the contrast-size dependence of magnetic features for var- 
ious positions on the solar disk can be fitted by parallel 
straight lines. A closer inspection of their Fig. 5 reveals 
that, especially for off-disk center features, the increase of 
contrast is sharper at the smallest areas, while no trend 
is observed at the largest are as. This is in agreement with 
Bovelet Wiehr I (|2003h and lHirzberger fc: Wiehr I (l2005h 



and also qualitatively with results obtaine d by analyses of 
full-disk images (see lErmolli et al. 20071 . and references 



therein). In order to investigate whether spatial resolu- 
tion of observations could explain such discrepancies, we 
studied the contrast-size relation of various clusters of flux 
tubes. Figure [TU shows the contrast-size relation for an 
isolated tube 70 km wide and for clusters of two, three 
and five tubes 70 km wide and 105 km apart at various 
positions on the solar disk and for spatial resolution 0.1" 
(top) and 0.3" (bottom). For vertical directions (/i = 1 
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Fig. 14. Size versus contrast of four clusters of flux tubes 
at various positions on the solar disk. Top: spatial resolu- 
tion 0.1". Bottom: spatial resolution 0.3". 

vary significantly with features size. At shallower lines of 
sight, the contrast increases with size. For /i = 0.8 and 
/i = 0.6 size dependence of contrast is larger for the small- 
est sizes while there is no dependence at the largest sizes. 
The increase of the contrast with size at inclined lines of 
sight is due to the reduced attenuation of radiation in the 
largest clusters. Moreover, since for each inclination, there 
is a finite number of flux tubes that radiation can cross 
(see §4), there is a size threshold beyond which contrast 
does not vary with size. This size threshold increases with 
the decrease of /i. An increase of contrast with size is ob- 
served even at vertical lines of sight with the reduction of 
spatial resolution, as illustrated in bottom panel of FiglMl 
This is due to the larger reduction in contrast value due to 
smoothing on the smallest features (see §4). Bottom plot 
shows that at each position on the solar disk, data can be 
well fitted with straight lines, as obtained in observational 



data bv Berger et all (|2007h . 



8. Summary and conclusions 

We have investigated thermal, photometric and geometric 
properties of various configurations of clusters of magnetic 
fiux tubes by two dimensional, static numerical simula- 
tions. The main results, obtained by the analyses of tubes 
with Wilson depression 60 km when in RE, are the follow- 
ing. 

At depths for which radiative heating is efficient, the 
temperature stratification within and around flux tubes 



tubes, but also of the tube number. In fact we found that 
the temperature inside the tubes increases with the num- 
ber of adjacent tubes in the optically thin part of the do- 
main and decreases at < log(r) < 1. At larger values of 
optical depth temperature increases again with the filling 
factor. The area around and between the tubes is cooled 
at 0.5 < r < 2, while is heated at greater depths. 

Since the temperature stratification determines photo- 
metric properties, contrast profiles of flux tubes also vary 
in clusters. At disk center contrast increases with the num- 
ber of flux tubes. Moreover regions in between the tubes 
are characterized by negative contrast, which is a signa- 
ture of the decrease of temperature, generated by the ra- 
diative cooling, in these areas. These dark features grad- 
ually disappear with the increase of the inclination of the 
line of sight. We found that for each cluster model there 
always exist an inclination beyond which clusters appear 
as compact monolithic features, as if they were isolated 
tubes. At these inclinations, the contrast of a cluster is 
larger than the contrast of a single flux tube element, and 
smaller than the contrast of an isolated tube with the same 
size and Wilson depression of the cluster. 

As for isolated flux tubes, off-disk center contrast pro- 
files of clusters are very asymmetric, with long tails toward 
the limb and dark lanes at disk center side. Nevertheless, 
projection onto the plane of the sky and the decrease of 
the spatial resolution (due for instance to both instrumen- 
tal and atmospheric degradation), largely makes profiles 
more symmetric. This finding is in agreement with mea- 
surements off-disk center of G -band bright points contrast 
profiles by Berger et alj (|2007l ) . These authors found quite 
symmetric profiles at all disk positions investigated. In 
both our models and observations, profiles at /i = 0.6 are 
more symmetric than profiles at /i = 0.4. This is true in 
the models for both isolated and clusters of fiux tubes 
(Fig. [13]). Reduction of spatial resolution also decreases 
the depth and size of dark rings and dark lanes, which 
disappear in the worst resolution cases investigated (0.3"). 
Moreover, dark lanes are not ubiquitous features of flux 
tubes. Since these are signatures of cooler regions within 
and around the tubes, they are very faint or absent in 
tubes models with_Jarger radiat ive heating; (see also mo d- 



Pizzo et al. I (|l993bh ;) 



els presented in ISteiner I (|20ut)l ) or 
or can be observed in some wavelengths a nd no t in oth- 
ers. This helps explain why Berger et al.l (|2007l ) did not 
find dark lanes in all the contrast profiles they observed 
off-disk center. 

CLVs of properties of magnetic bright features are 
also affected by clustering and reduction of resolution. 
The contrast-CLV is for instance steeper for isolated flux 
tubes than it is for clusters of tubes. Since most of feature 
identification techniques on images are based on intensity 
threshold, isolated small (and therefore less brilliant) fea- 
tures are less likely to be detected. Measurements at the 
limb are therefore biased toward larger tubes and clusters 
and the resulting CLV is flatter than the one obtained 
for single structures (see Fig. [7j). For the same reasons. 
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[TQ|) . This would explain CLV shapes obtain ed by recent 



obser vations at the extreme limb (Hirzbe rger fc Wiehr 



2QQ5h . Reduction of resolution also shifts the peak of the 



contrast-CLV toward lower values of /i, especially for tubes 
models with marked dark lanes, and makes the CLVs flat- 
ter (see Fig. [8j). In this context, we notice that Berger 
et al. 2007, who analyzed better resolution images with 
respect to other authors, found the peak of the CLV of 
contrast to occur closer to the disk center with respect to 
some previous works (see Berger et al. 2007 and references 
therein). 

Finally, we have investigated the size-contrast relation 
for clusters of various numbers and spatial resolutions. We 
found a slight increase of contrast with size at disk cen- 
ter, which has to be partially ascribed to the increase of 
temperature with the number of adjacent tubes and par- 
tially to the larger effects of smoothing on smaller features. 
Increase of contrast with size (and therefore with the num- 
ber of flux tubes) is found for inclined lines of sight due to 
the reduced opacity inside the tubes. Since for a given in- 
clination of the line of sight only a finite number of tubes 
is crossed by radiation (see §4), there is a size thresh- 
old beyond which contrast does not vary. This threshold 
increases with the inclination of the line of sight. With 
the decrease of spatial resolution, even for vertical lines of 
sight the contrast increases with size, the threshold size 
value increases, and the contrast-size curves obtained at 
various lines of sight can be fitted by straight parallel lines. 
These findings are in agreement with recent contrast-siz e 
measurement s obtained on both h igh ([Berger et al.ll20Q7l ) 
and medium (lErmolh e t"ani2007h resolution images. 

Results presented in this work suggest that measure- 
ments of properties of magnetic bright features are sig- 
nificantly influenced by spatial resolution, image analysis 
techniques, and the wavelength of observation. Geometric 
properties of contrast profiles, such as size and asymme- 
tries, are particularly affected. These findings can par- 
tially explain the discrepancies presented in the literature. 
Moreover, some observational results can be better inter- 
preted by considering observed bright magnetic features as 
aggregation of smaller elements, rather than a monolithic 
entity. 

The results presented in this work should be properly 
taken into account when trying to derive physical proper- 
ties of magnetic flux tubes by observations. 
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